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Light and Strange Quark Masses with Dynamical Wilson Fermions. 

Henning Hoeber (for the SESAM collaboration*) 8, 

a Department of Physics, University of Wuppertal, Gaussstr. 20, 42119 Wuppertal, Germany; 
HLRZ c/o Jiilich Research Centre and DESY, Hamburg, 52425 Jiilich, Germany. 

We determine the masses of the light and the strange quarks in the MS-scheme using our high-statistics lattice 
simulation of QCD with dynamical Wilson fermions. For each of our three sea quarks we have analyzed our data 
at five different values of the valence quark mass, enabling us to parameterize our fit results in the (m sea , m valcncc ) 
plane. For the light quark mass we find milS^(2 GeV) = 2.7(2) MeV, which is lower than in quenched simulations. 



Applying a new method, which we propose to extract the strange quark mass in a sea of two dynamical light 
quarks, we obtain mf^ nso (2 GeV) = 140(20) MeV. 



1. Introduction 

More than 25 years after the introduction of 
quarks as the constituents of hadrons the absolute 
values of the light quark masses, as given in the 
current Particle Data Book (!]] , are still among the 
most poorly known fundamental parameters of 
the Standard Model. Chiral perturbation theory 
(xPT) has been very successful in fixing the ratios 
of light quark masses || but the overall scale is 
beyond its possibilities. QCD sum-rules || have 
produced a wealth of predictions but in a large 
range of values and systematic errors are hard 
to pin down. Given the shortcomings and prob- 
lems of these two methods, the lattice, as a first- 
principles non-perturbative method, would ap- 
pear to be the most promising candidate for a reli- 
able extraction of the absolute quark mass values. 
In this talk I present results of the SESAM collab- 
oration for the masses of the light and strange 
quarks obtained from the currently largest and 
statistically most precise lattice QCD simulation 
with 2 flavours of dynamical Wilson fermions j| . 

2. Method and Status 

Calculating quark masses on the lattice con- 
sists of three steps: firstly, a hadron mass is 
matched to its experimental value to fix the bare 
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lattice quark mass under consideration. Typi- 
cally, the pseudoscalar meson mass, having the 
smallest statistical lattice-error of all hadrons, is 
used. This step also invokes xPT, often only to 
lowest order, to relate the pseudoscalar mass to 
the bare quark mass. Secondly, the lattice scale 
is set in a similar manner by matching some ex- 
perimentally well known quantity such as the p 
mass or the 15 — 2S energy level splitting in the 
T spectrum. Finally, perturbation theory is used 
to convert the bare quark mass to that of the M S 
scheme. 

Quite a large number of simulations have been 
carried out in this manner in the quenched ap- 
proximation of QCD. Internal fcrmion loops are 
totally neglected in these simulations for want of 
sufficient compute power. Recent compilations 
and interpretations of these results are given in 
|U and §. The authors of §,§ have each 
attempted to produce quark mass values in the 
limit of vanishing lattice spacing. Although the 
analysis differ greatly in approach and method 
the overall result is very much the same: the 
isospin symmetric light quark and the strange 
quark masses are both predicted by quenched lat- 
tice QCD to be smaller than the sumrule results 
§. For full QCD, very few results have been ob- 
tained so far but these would seem to indicate 
even smaller quark masses; however, no reliable 
continuum extrapolation can be performed so far. 
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3. SESAM analysis 



SESAM has recently completed the generation 
of gauge configurations with two dynamical Wil- 
son fermions at three different values of the dy- 
namical mass and for a lattice extent of 16 3 x 32. 
The bare coupling of ft = 5.6 corresponds to a 
cutoff a" 1 ~ 2.33(6)GeV (matching the p)|. 200 
lattices are analysed per sea-quark mass, giving 
m„/m p = 0.841(6), 0.755(8), and 0.69(2). Cor- 
relators are evaluated for 5 values of the valence 
quark mass at each sea quark mass value. 

3.1. Light Quark Mass 

To extract the light quark mass from the pseu- 
doscalar (PS) and vector (V) meson mass ratio 
m^/mp we attempt to fit the masses of mesons 
with valence quark masses equal to sea quark 
masses (characterized by the argument ss) ac- 
cording to first order xPT: 



'PS 



(ss) 



mv{ss) = m crlt +bm 2 



PS 



•(1) 



Figure [l] shows that this ansatz works very well 
(although a quadratic term in the vector data 
cannot be excluded). The condition 



mps(m l q ) _ 
my(m^) \m 



= 0.1785 



(2) 



cxp 



yields 



MS 



(2GeV) = -(m„ 



m d ) 



2.7(2) MeV . (3) 



The lattice result has been converted to the MS 
scheme using boosted 1-loop perturbation theory 
H and the scale was set by the p-mass. 

3.2. Strange Quark Mass 

So far, we have used our data with equal sea 
and valence quark masses to extract the light 
quark mass in a sea of light quarks. To simu- 
late the strange quark we need to introduce va- 
lence quarks that are unequal to the dynamical 
quarks. The meson masses with strange valence 
quarks can then be extrapolated to the sea of light 
quarks. To this end, encouraged by eq. 0, we 



2 If, instead, we match the 15 — 25 splitting in the T, the 
scale turns out to be ~ 2.48(14)GeV. 
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Figure 1. m| s ss (left) as a function of and 
mv,ss as a function of m sea (in lattice units). 
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make the all-linear ansatz 
m 2 PS (sv) 



l PS 



(vv) 



c m n 



c"m s q ea 



d'm v q al , 
-d"(rrC 11 



(4) 
.(5) 



The arguments sv and vv denote the valence 
quark content so that, for example, vv means nei- 
ther of the valence quark masses is equal to the 
sea quark mass. 

Equations [| and || seem to show that four ad- 
ditional parameters are needed to parametrize 
the pseudoscalar meson mass in the entire (m sea , 
m valenco ) plane. However, the consistency of 
equations |^ and |^ with eq. [I] can be used to reduce 
the number of unknown parameters to 3 and sim- 
ilar arguments hold for the vector particle. Thus, 
we attempt fits to the combined parametrization 



TO PS,5S 

2 

PS.sv 



m 

,2 




val . 



'PS.vvM 

and similarily for the vector meson. Figure | dis- 
plays our data and the combined fits according 
to this ansatz. Again, we find the parametriza- 
tion to work very well but cannot exclude small 
contributions from higher order x?T. 

The strange quark is now determined from the 
(redundant) system of conditions 



m K , 



= 1.16, 



exp 



( UUL] 

\ m p ) t 



= 0.642 



m 



1.326 , 



which has to be evaluated at m q ea — m l q . All 
three conditions lead to consistent results for m s ; 
we average and include the spread into the statis- 
tical error. After conversion into the MS scheme 
we obtain 



140(20) MeV . 



(7) 



4. Discussion 



Our dynamical results can be compared to the 
analogous quenched values at similar lattice spac- 
ing (corresponding to /3 qU onchcd = 6.0): 



mg*(2GeV) = 
^ ngc (2 GeV) 

MS v ; 



5.5(5) MeV 
= 166(15) MeV 
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Figure 2. Simultaneous fit of all pseudoscalar and 

vector data to eq. [| Symbols: * = m ss ; 

O = m sv — ■ — ; □ = m m . All plots in lattice 

units. 
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Errors due to the finite lattice spacing and the 
finite volume are not included in (|L and (||). 

Compared to the quenched result (at the same 
cutoff) , which is in good agreement with previous 
lattice calculations j(| , we observe a much smaller 
dynamical light quark mass, whereas the strange 
masses are compatible within errors. It is impor- 
tant to keep in mind that the dynamical light and 
strange quark masses have been obtained with 
two slightly different methods: the light quark 
is obtained in two-flavour QCD whereas evaluat- 
ing the strange quark mass necessitates using va- 
lence quarks unequal to the sea quarks. This may 
explain why the ratio of light to strange quark 
masses, so tightly constrained by X-PT, is not re- 
produced by the two dynamical values of eqs. [j| 
and @. 

To illucidate this matter further, it is instruc- 
tive to examine the values of the light quark 
masses at fixed sea quark masses. The surpris- 
ing result is that extrapolating the valence quark 
to is not sufficient to make the pseudoscalar 
masses mp s (sv) and m 2 PS (vv) vanish at fixed 
sea quark masses. Instead, we need to tune the 
valence quark masses to negative values. The 
reason for this can be traced back to the defi- 
nition of the quark masses, which, for Wilson- 
like theories, where chiral symmetry is broken, 
are definded relative to the critical value of the 
hopping parameter k c . In the dynamical theory, 
the critical point is taken at vanishing sea quark, 

Kc = K c (m sc& — 0), m q = i f~ — -l-J , where the 
mass of the pseudoscalar truely vanishes. If in- 
stead, we define a variable, sea quark mass depen- 
dent critical kappa at each sea quark value we find 
much larger quark masses (5.7(4), 5.6(3), 5.4(3) 
MeV) which are similar to the quenched quark 
mass. The situation is displayed schematically in 
figure ||. Ai is the value of the light quark mass 
we obtain in the full theory from m| s ss whereas 
values such as A2 are obtained from fixed sea- 
quark measurements (partially quenched). We 
see from figure || that the light quark masses from 
finite sea-quark mass values do not extrapolate 
to the value Ai. Thus, to extract the light quark 
mass, it it necessary to work with rrip S 5S ; the data 
and m| s vv can be used as a consistency 
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Figure 3. Schematic plot of light quark mass de- 
terminations (see text). 

check but are not sufficient to obtain m hght . 

To consolidate our results we are presently cal- 
culating the quark masses using an alternative 
method in which the bare quark mass is defined 
by means of the PCAC relation, 

mq ~ 2(P(x)P(0)) ' [> 

where A is the axial current and P the pseu- 
doscalar density. This method may be less sen- 
sitive to renormalisation effects as the previous 
one, thus giving an important cross-check. 
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